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ABSTRACT 

We explore the possibility of explaining the properties of the Be/X-ray pulsars observed in 
the Small Magellanic Cloud (SMC) within the magnetic levitation accretion scenario. This 
implies that their X-ray emission is powered by a wind-fed accretion on to a neutron star (NS) 
which captures matter from a magnetized stellar wind. The NS in this case is accreting matter 
from a non-Keplerian magnetically levitating disc which is surrounding its magnetosphere. 
This allows us to explain the observed periods of the pulsars in terms of spin equilibrium 
without the need of invoking dipole magnetic fields outside the usual range ~ 10" - lO" G 
inferred from cyclotron features of Galactic high-mass X-ray binaries. We find that the equi¬ 
librium period of a NS, under certain conditions, depends strongly on the magnetization of 
the stellar wind of its massive companion and, correspondingly, on the magnetic field of the 
massive companion itself. This may help to explain why similar NSs in binaries with similar 
properties rotate with different periods yielding a large scatter of periods of the accretion- 
powered pulsar observed in SMC and our galaxy. 

Key words: accretion, accretion discs - pulsars: general - stars: winds, outflows - X-rays: 
binaries. 


1 INTRODUCTION 

The Small Magellanic Cloud (SMC) contains a large number of 
accreting pulsars in high mass X-ray binaries (HMXBs). These 
sources, being at a well known, virtually uniform distance and with 
small interstellar absorption, constitute an ideal sample for popula¬ 
tion studies of neutron star (NS) binaries (Haberl & Pietsch 2004; 
Shtykovskiy & Gilfanov 2005; Laycock et al. 2010). 

Klus et al. (2014) have recently reported parameters of 42 
Be/X-ray binaries in the SMC observed with the Rossi X-ray Tim¬ 
ing Explorer (RXTE) satellite over a time span of 14yr. These sys¬ 
tems, most of which are transients, contain pulsars in which the 
X-ray emission is powered by wind-fed accretion on to a mag¬ 
netized NS. The average spin periods of these pulsars are in the 
range ~ 2.37 - 1323 s and change at an average rate \P\ ~ 
(0.02 - 620) X 10“^ s yr“'. Their orbital periods range from about 4 
to 500 d. 

Klus et al. (2014) estimated the magnetic fields of these NSs 
considering the situation in which they accrete matter from a Ke- 
plerian disc and rotate close to their equilibrium spin period. In 
this case, the estimated fields of the majority of the stars (including 
all those with > 100 s) turn out to be over the quantum criti- 
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cal value = nP^c^lefi == 4.4 x 10*'^ G. Here and e are the 
mass and electric charge of an electron, c is the speed of light and 
Ti is the reduced Planck constant. Alternatively, if these NSs are not 
close to spin equilibrium, their inferred magnetic fields are smaller 
than ~ 10*** G. Such results are rather unexpected since the ma¬ 
jority of NSs in our Galaxy, including accretion-powered pulsars 
in HMXBs in which the magnetic field is measured from the cy¬ 
clotron resonance scattered features, have magnetic fields between 
10" and 10" G (see, e.g., Revnivtsev & Mereghetti 2015). An at¬ 
tempt to invoke currently used quasi-spherical accretion scenarios 
did not help much to improve the situation either, leading to esti¬ 
mated magnetic fields in excess of lO" G. 

It is not unusual that the magnetic fields evaluated from the 
spin parameters of NSs in Galactic HMXBs significantly exceed 
those measured through observations of their cyclotron resonance 
features. As recently indicated by Ikhsanov et al. (2014), this in¬ 
consistency may reflect an oversimplification of currently used 
wind-fed accretion scenarios, in which the magnetic field of the 
matter captured by the NS from its environment is neglected. The 
incorporation of the fossil magnetic field of the accreting matter 
into the model leads, under certain conditions, to a different accre¬ 
tion regime which is referred to as magnetic levitation accretion. In 
this scenario, which is briefly outlined in the next section, accre¬ 
tion occurs through a non-Keplerian magnetically levitating disc 
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(ML-disc) and the maximum possible torque exerted on the NS 
significantly exceeds that previously evaluated in traditional non¬ 
magnetic scenarios. This leads to a new expression for the equilib¬ 
rium period (see Sections), which allows us to explain the observed 
values of the spin and orbital periods of the SMC pulsars with¬ 
out invoking magnetic fields outside the canonical range of 10*'- 
10*3 G (Section 4). 


2 MAGNETIC LEVITATION ACCRETION 

We consider an HMXB, with orbital period composed of a 
magnetized NS rotating with spin period and a massive early- 
type star, which underfills its Roche lobe and loses matter through 
a stellar wind. The X-ray emission of the system is powered by 
wind-fed accretion on to the NS. This implies that the NS cap¬ 
tures matter from the wind at a rate 911 ^ 9jl(. = where 

= 2GMns/i'3^j is the Bondi radius, is the NS mass, v^i is its 
velocity in the frame of the wind, and po = pfr^,) is the density of 
matter in the region of interaction. The captured matter moves to¬ 
wards the NS forming an accretion flow, which interacts with the 
stellar magnetic field and confines the magnetosphere within the ra¬ 
dius r^, where the flow pressure is balanced by the pressure of the 
NS magnetic field. The accreting matter penetrates into the field 
at the magnetospheric boundary and, finally, falls on to the stellar 
surface at the magnetic pole regions by moving along the magnetic 
field lines. 

The analysis of this accretion scenario reported by 
Ikhsanov & Finger (2012) indicates that the structure of the ac¬ 
cretion flow, as well as the appearance of the accretion-powered 
source, strongly depend on the physical conditions of the matter 
captured at the Bondi radius, which, in general, possesses some an¬ 
gular momentum and magnetic field. Under these conditions, the 
structure of the accretion flow beyond the magnetosphere can be 
treated in the following basic approximations: (i) a spherically sym¬ 
metrical or quasi-spherical flow, (ii) a Keplerian disc, and (iii) an 
ML-disc. A key parameter which allows us to determine which of 
these situations applies is the relative velocity of the NS with re¬ 
spect to surrounding matter, Vrei. 

The spherically symmetrical accretion occurs if the gas sur¬ 
rounding the NS does not possess angular momentum and is non- 
magnetized. The captured matter in this case moves towards the 
NS in the radial direction with the free-fall velocity, Vff(r) = 
(2GM„s/r)*^3 ijj a spherically symmetrical fashion. Its density 
scales with the radius as p(r) = 9)1/ (4tTr^Vff'j and the ram pressure 
is 

£ram('') = P(f)Vff(r) oc ^3/2^ (j) 

The minimum distance to which the spherical flow can approach a 
NS with dipole magnetic moment pis r > r^ (Arons & Lea 1976; 
Eisner & Lamb 1977), where 

- 

\9JI(2GM„,)'/2 

is the Alfven radius which is defined by equating the ram pressure 
of the free-falling gas, with the magnetic pressure due to dipole 
magnetic field of the NS, Pm = p^! (27Tr^'j oc 

The structure of a non-magnetized accretion flow deviates 
from spherically symmetrical if the captured matter possesses an¬ 
gular momentum. The process of mass accretion in this case is ac¬ 
companied by the accretion of angular momentum, which in a bi¬ 
nary system with the orbital angular velocity occurs 



at the rate J = (Illarionov & Sunyaev 1975). Here ^ is 

the parameter accounting for dissipation of angular momentum due 
to density and velocity gradients in the accreting non-magnetized 
gas (see e.g. Ruffert 1999, and references therein). The angular 
velocity of matter in this so called quasi-spherical accretion flow 
scales with the radius as 

nr(r) = fn„,b (y)\ (3) 

and the azimuthal component of the dynamical pressure of the flow, 

£.flr)=pof- Q.l^, '•q (y) ocr-''/^, (4) 


increases more rapidly than the ram pressure. The accretion pro¬ 
ceeds in a quasi-spherical fashion up to the circularization radius 


Girc — 



(5) 


at which condition E-flr) = £ram('') is satisfied. The angular veloc¬ 
ity of the accreting matter, Gf(r), at this radius reaches the Keple- 
rian angular velocity, = I' /2GAf„sj , and the accretion flow 
switches into a Keplerian accretion disc (Pringle & Rees 1972; 
Shakura & Sunyaev 1973). 

The accretion picture may differ from that presented above if 
the matter captured by the NS at the Bondi radius is magnetized. As 
long as the Alfven velocity, Va(r) = 6f(r)/ \Anp{r)\^^^, in the accret¬ 
ing matter is smaller than the free-fall velocity, the magnetic field 
of the accreting gas, Bf, does not influence significantly the flow 
structure. Therefore, the captured matter initially follows ballistic 
trajectories forming a quasi-spherical accretion flow in which the 
angular momentum and magnetic flux are almost conserved. How¬ 
ever, the magnetic pressure in the free-falling gas increases rapidly, 

c . ^ -4 c.c 

£m(f) = -r— oc r , (6) 

8;r 


and reaches the ram pressure at the Shvartsman radius (Shvartsman 
1971) 





(7) 


which is defined by equating the magnetic to ram pressure. Here 
flo = j8(r(,) = Snpocl^l is the ratio of thermal to magnetic pres¬ 
sure in the captured matter at the Bondi radius, Cso = Csfr^,) is the 
sound speed and Sfo = Bflr^) is the magnetic field in the accreting 
matter. The flow at this radius is converted into a non-Keplerian 
ML-disc in which the accreting matter is confined by its own 
magnetic field (for discussion see Bisnovatyi-Kogan & Ruzmaikin 
1976; Igumenshchev, Narayan, & Abramowicz 2003). The accre¬ 
tion inside the ML-disc proceeds on the time-scale of the magnetic 
flux dissipation and, in the general case, can be treated in a diffu¬ 
sion approximation. The value of the effective diffusion coefficient 
strongly depends on the configuration of the magnetic field in the 
disc and its stability. It ranges from the Bohm diffusion coefficient, 
if the magnetic field annihilation in the disc is governed by dissipa¬ 
tive instabilities and magnetic reconnections, up to much larger val¬ 
ues if the field configuration is interchange unstable (for discussion 
see Igumenshchev 2009; Tchekhovskoy, Narayan and McKinney 
2011 ; Dexter et al. 2014, and references therein). 

Thus, a quasi-spherical accretion on to a NS in a wind-fed 
HMXB occurs if the Alfven radius is larger than both the circu¬ 
larization and the Shvartsman radii. The angular momentum and 
magnetic field of the accreting material in this case are too small to 
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significantly influence the flow structure before the ballistic trajec¬ 
tories of the accreting matter are truncated by the magnetic field of 
the NS. Solving the inequality > maxjrcirc,f?shl for Vrei one finds 

I Vcn for rcirc, 

( 8 ) 

for ^ ^circ) 

where (Ikhsanov 2007) 

v„ - 160km (9) 

and (Ikhsanov & Beskrovnaya 2012; Ikhsanov & Finger 2012) 

v„a-460kms-‘x/3„'/V3o'^^"‘^'^^™r 

Here /J 30 = /r/ (lO'^® Gcm^) is the dipole magnetic moment (// = 
(l/2)6nsl^ns) of a NS with surface magnetic field B„s and radius 
m = M„,J1.4Mo, 'Jit ,5 = W/(l 0 >'gs-‘), F 50 = F„b/50d, 
and the parameter ^ 0.2 = ^/0.2 is normalized according to Rutfert 
(1999). The radius of the magnetosphere, within this scenario, is 
comparable or slightly exceeds the Alfven radius (Arons & Lea 
1976; Eisner & Lamb 1977), and accretion beyond the magneto- 
spheric boundary proceeds in the form of a free-falling gas or a 
hot turbulent spherical envelope (Lamb et at. 1977; Ikhsanov 2001, 
2003). 

A Keplerian accretion disc in a wind-fed HMXB can form if 
the circularization radius exceeds both the Alfven and Shvartsman, 
radii. Solving inequality r^irc > maxjr^,Rsi,l for v,.ei one finds 

( Vcr, for R* < 

( 11 ) 

[ Vca, for R,h > 


(for discussion see Ikhsanov 2012; Ikhsanov et at. 2014). Here 
is the proton mass, kg is the Boltzmann constant and To is the tem¬ 
perature of the matter at the inner radius of the disc. The parameter 
Og = D^glDg is the ratio of the effective coefficient of diffusion 
of the accreting matter into the stellar field at the magnetospheric 
boundary, D^g, to the Bohm diffusion coefficient, which in the con¬ 
sidered case can be expressed as Dg = ckg Tgr^.^H'ileij). Numerical 
simulations of the diffusion process and measurements of the rate 
at which the solar wind penetrates into the magnetosphere of the 
Earth suggest that ~ 0.01-1 (Gosling et at. 1991). The matter 
being diffused into the stellar field flows along the magnetospheric 
field lines and reaches the NS surface at the magnetic pole regions. 

The possibility of explaining the parameters of the SMC 
Be/X-ray pulsars within the quasi-spherical and the Keplerian 
disc accretion scenarios has been already discussed by Klus et al. 
(2014). Here we explore the possibility that these pulsars accrete 
matter from an ML-disc. In the next section we evaluate the equi¬ 
librium period of an NS in this scenario and present the expected 
Corbet (P^ versus Porb) diagram in Section 4. 


3 EQUILIBRIUM PERIOD 

The equation governing the spin evolution of an NS accreting mat¬ 
ter from an ML-disc reads 

iTtlv = K.^ + Kg + Ke, (17) 

where I is the moment of inertia and v = l/P^ = tUs/2;r is the spin 
frequency of the NS. The first term in the right-hand side of this 
equation, 

K, = me(r^,)=^muj,ri„ (18) 


where (Ikhsanov, Kim & Beskrovnaya 2015) 

Vca - (12) 

Finally, a non-Keplerian ML-disc can form in a wind-fed 
HMXB if the Shvartsman radius is larger than both, the Alfven and 
circularization radii. Solving inequality > maxjr^, rcird for Vrei 
yields 

Vca Ael < Vma. (13) 

The inequality Vca < Vma is satisfied if /3o < /Imax, where 


is the rate at which angular momentum is transferred to the NS 
by the matter flowing inside the magnetosphere. Here 
rma X (rmaWs) Is the specific angular momentum of matter in the 
magnetopause at the rotational equator. This matter, which pene¬ 
trates into the stellar field at the magnetospheric boundary, coro¬ 
tates with the NS and flows towards its surface along the field lines. 
Ka has a positive sign and represents the minimum possible spin-up 
torque exerted on to a star which is surrounded by a magnetosphere 
of radius r,na and accretes matter on to its surface at a rate Sli. 

The second term (Ikhsanov 2012; Ikhsanov et al. 2014), 


ySmax - 164 m 


F 30 


i 1 

2 ^ hf 1 



1 ) 

\ 10km s “'} 


(19) 


This indicates that the magnetic levitation accretion scenario oc¬ 
curs if the fossil magnetic field in the matter captured by the NS is 
6f(G) > where 


( 2mvLci 


( 




1/2 


^ 6 X 10 G X 


(15) 


/ Vrei 

VlOOkms^*/ 


5 ° ‘5 (lookms- 


The matter in an ML-disc approaches the NS until the radius 

f cm] frr^/'y/‘3(GM„G‘/‘^ 
ll6V2.fcBi ’ 


accounts for the angular momentum exchanged between the star 
and the disc at the magnetospheric boundary (here ki is a dimen¬ 
sionless parameter of order unity). The sign of this term depends 
on the ratio between the angular velocity of matter at the inner ra¬ 
dius of the disc, Qf(r^ii), and the angular velocity of the NS. 

The last term, accounts for the angular momentum ex¬ 
changed between the NS and the disc at radii r > r^a. Under the 
conditions of interest (dGf/dr < 0), this term is smaller than Kg 
at least by a factor ~ (rjr^^)^^'^. Therefore, its contribution is rela¬ 
tively small and we neglect it in the following. 

The equilibrium period, Peq, is defined by equating the total 
torque exerted on the NS (right-hand side of equation 17) to zero, 
which, using equations (18) and (19), yields 


where the pressure exerted by the disc on to the stellar magneto¬ 
sphere is equal to the magnetic pressure due to the NS dipole field 
and the rate of diffusion of the accreting matter into the stellar field 
is equal to the mass capture rate by the star from its environment 


— ^f(^ma) 


1 - 





where Rf(rmJ = 27r/Gf(r,„a). 


( 20 ) 
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This equation suggests that the angular velocity of the NS in 
spin equilibrium exceeds the angular velocity of the matter at the 

inner radius of the disc by a factor of |l - (1/ Vz^t) (''ma/'li,) ^ |- 
The torque K\, in this case has a negative sign and tends to spin- 
down the star. It, however, is compensated by the spin-up torque K.^ 
and thus, the total toque exerted on the NS is zero. 

Evaluating the ratio 

^ 10-3 tt;/13 ^-5/13 ,„10/13 3il-l/13y^-7/13^ ^21) 

one finds that the expression in the brackets in the conditions of 
interest (10^3 <k < 1) is close to unity (here = Tq/IO^K is 
normalized to the temperature of optically thick gas irradiated by 
the pulsar emission, see, e.g. Hickox, Narayan and Kallman 2004). 
This indicates that the angular velocity of the NS in spin equilib¬ 
rium is close to the angular velocity of the matter at the inner radius 
of the disc. This allows us for simplicity to approximate the equi¬ 
librium period of the NS by 

The analysis of the transport of angular momentum across the 
ML-disc is rather complicated and beyond the scope of this paper. 
The angular momentum extracted from the NS can be stored in the 
disc (Sunyaev & Shakura 1977; D’Angelo & Spruit 2012) or/and 
transferred out from the disc by Alfven waves. Here we simply 
assume that the angular velocity of the matter in the disc scales 
with the radius as flf(r) oc (RshlrY, where y is a free parameter, 
which in the general case is limited to* y < 2. The angular velocity 
of matter at the inner radius of the disc can be expressed as 

Df(rma) = nf(7?*)(^) , (22) 

where 

a,(l?sh)-^f2orb|^) (23) 

is the angular velocity of matter at the Shvartsman radius. This 
is derived by taking into account that the magnetic field of the 
accreting matter does not significantly influence the flow struc¬ 
ture beyond the Shvartsman radius and the matter in the region 
Psh < r ^ moves towards the NS quasi-spherically. 

Putting equations (7), (16) and (23) into equation (22) and tak¬ 
ing into account that the relative velocity in the considered case is 
''rei = ((Vma, with ^ < 1, One finds 

/I 3107/3\ 

Pf(rma) = 0.5 (24) 

where 

( cin^OL 

A,, = (25) 

and 


fli = (208 -50y)/455, 
fl 2 = (lOy-H 104)/455, 
03 = (10y-8)/3. 


' The value y=2 corresponds to conservation of angular' momentum, y= 1.5 
to Keplerian rotation, y = 1 to rotation with constant linear velocity, and y=0 
to solid body rotation. 


04 = (100y-416)/455. 


Finally, the value of the dimensionless parameter (in the gen¬ 
eral case lies in the interval < ^ < 1, where 


^min — 


^ns/^ma> 

for Vns > Vca 

^ca/^^ma> 

for Vns < Vca 


and t’ns is the orbital velocity of the NS. 


(26) 


4 APPLICATION TO THE SMC PULSARS 

In this section we explore the possibility of explaining the spin 
properties of the SMC Be/X-ray pulsars in terms of the mag¬ 
netic levitation accretion scenario described above. We assume that 
the NSs in these systems rotate close to the equilibrium period 
Peq - Pf(rma) giveu by equation (24) and accrete matter from 
an ML-disc. This implies that the relative velocity of the NSs in 
the frame of stellar wind of their companions meets the condition 
Vcr < Vrei < v^a and/ffl < /6max expressed by equation (14). Since the 
accretion process in an ML-disc is fully controlled by the magnetic 
field of the accreting matter itself, we consider the case of rigid 
rotation of the disc by setting y = 0. 

The minimum possible value of the equilibrium period, P^m, 
which an NS can acquire within the magnetic levitation accretion 
scenario can be derived from equation (24), by setting fio = Anax 
and, correspondingly, ^ = 1. This yields 

(27) 

where DJlis is the average mass accretion rate normalized to 
10*3 gs-'. 

In Fig. 1 the values of are compared to the observed val¬ 
ues of Ps and Porb for representative values of /j (103** G cm3 
lines; lO^* G cm3 p^es) and 9il (10*“* and 10*^ g s“*). The 

corresponding lines indicates the minimum possible equilibrium 
period of an NS which captures matter from a relatively slow, mod¬ 
erately magnetized stellar wind in the magnetic levitation accretion 
scenario. 

The maximum possible equilibrium period, Pmax, can be eval¬ 
uated from equation (24) by setting /Bo = 1 and ^ = ^^in, which for 
the case v„s < Vca is (see equation 26): 

- 0.17/li^/35 ^6/35 3jj-3/35 ^3/7 ,„3/35 ^^24/35 ^-3/7^ (28) 

Putting this value into equation (24) one finds 

PZ. = 20 s X ^ ^'*^/’ cf(29) 

The function PZx = 7’max(7’orb) is shown by the dot-dashed line in 
Fig. 1. The maximum equilibrium period of an NS which accretes 
from a magnetized slow wind does not depend on the magnetic field 
of the NS itself and on the X-ray luminosity of the pulsar. However, 
it is a strong function of the sound speed in the surrounding matter 
and, therefore, it can exceed the value given by equation (29) if the 
temperature of the stellar wind exceeds 10“* K. 

Fig. 1 shows that practically all the pulsars reported by 
Klus et al. (2014) have spin-period values smaller than P^^^ = 
Pmax(Porb) and lying above the lines of P^in corresponding to 
/j ~ 103^-3** G cm3 accretion rates consistent with the long¬ 
term average luminosities of these sources. This suggests that the 
observed spin periods can be explained within the magnetic levita¬ 
tion accretion scenario with surface magnetic fields of the NSs in 
the canonical interval ~ 10*'-10*3 G. 
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Figure 1. Corbet diagram of the Be/NS pulsars in the SMC, with lines of maximum and minimum equilibrium period in the ML-disc accretion hypothesis. 
The solid lines are the minimum equilibrium period for fi = 10'^® G cm^ and 9ji = lO'^ and lO'^ g s“'. The dashed lines are the minimum equilibrium period 
for = 10^® G cm^ and 'Jji = 10*'* and 10*^ g s“*. The dash-dotted line indicates the maximum equilibrium period (independent of fi and 3)i). 


5 DISCUSSION 

We find that the observed spin periods of the SMC Be/X-ray pulsars 
are in a range consistent with the values expected for the equilib¬ 
rium periods of NSs with magnetic fields of 6 ~ 10**“'^ G and 
accreting from an ML-disk. The assumption that the NSs in these 
systems are spinning close to an equilibrium value is supported by 
the fact that they show alternate episodes of spin-up and spin-down, 
which do not change significantly on the long term. We believe 
that the equilibrium period is set by the average accretion rate ex¬ 
perienced by these NS during the long time intervals of quiescence 
(or low X-ray luminosity) between their bright outbursts. In fact 
these sources spend most of the time in such low-'Jj! conditions, re¬ 
sulting in X-ray luminosities well below the RXTE sensitivity limit 
of ~ 10^^ ergs“‘. For this reason we have adopted in Fig. 1 values 
of 3jt corresponding to luminosities of ~ 10**'*“^^ erg s“*. 

The equilibrium period of an NS which accretes from a ML- 
disc tends to increase with the orbital period of a binary system 
(see equation 24). This, in particular, can be a reason for a lack of 
pulsars in the lower-right part of the Corbet plot. On the other hand, 
Fcq also depends on several other parameters, such as the magnetic 
field and relative velocity of an NS, the mass accretion rate, and 
physical conditions in the stellar wind with which the NS interacts. 
The great diversity of possible combinations of these parameters 
is responsible for the large scatter of the observed spin periods in 
the Corbet plot. In this accretion regime, the magnetic field of the 
massive star plays an important role, with stronger fields leading to 
longer equilibrium periods of the pulsar. 

Spectropolarimetric observations of O/B-type stars give evi¬ 


dence for a relatively strong magnetization of these objects (see, 
e.g. Walder, Folini & Meynet 2012, and references therein). The 
strength of the large-scale field at the surface of several of these 
objects has been measured in the range ~ 500 - 5000 G, and in 
some cases beyond 10 kG (Hubrig et al. 2006; Martins et al. 2010; 
Oksala et al. 2010). Some of the early type stars are surrounded by 
X-ray coronae which indicate the magnetic activity of these ob¬ 
jects (Schulz et al. 2003). As most of these stars rotate relatively 
fast (see e.g. Rosen, Krumholz & Ramirez-Ruiz 2012, and refer¬ 
ences therein), the magnetic field in the wind is dominated by the 
toroidal component which scales with the radius oc r“*. 

A similar situation is realized in the solar wind in which the 
magnetic field at a distance of 1 au is 6sw ~ 10“^ G and the param¬ 
eter ji is close to the equipartition value, p ~ \ (Mullan & Smith 
2006). Following this similarity one can suggest that the surface 
large-scale magnetic field of massive stars in the considered sys¬ 
tems is a factor of ~ 50 larger than the surface large-scale 

magnetic field of the Sun and can be as large as a few hundred 
Gauss. A smaller magnetization of these stars cannot be also ex¬ 
cluded if the dynamo action applies in their outflowing discs. The 
fact that periods of the considered pulsars are much shorter than the 
maximum possible period predicted by our model may indicate that 
the average value of Pq exceeds unity and hence, the stellar wind of 
early spectral type stars is less magnetized than the solar wind. The 
observed range of periods of the pulsars (~ 1 - 1000 s) in this case 
can be explained in terms of variation of Pq parameter from system 
to system within an order of magnitude. 

Our study confirms the conclusion of Klus et al. (2014) that all 
of the considered pulsars are situated in a relatively slow wind. This 
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is consistent with current views on the mass outflow process of Be 
stars, in which the stellar wind at the equatorial plane is dominated 
by a dense outflowing disc. The radial velocity of matter in the disc 
is comparable or even smaller than the orbital velocity of the NS 
(Okazaki & Negueruela 2001). 

We finally note that our results apply also to the Be/X-ray pul¬ 
sars in our Galaxy, which show a distribution in the Poa-Ps similar 
to that of the SMC sources. This supports the view that the proper¬ 
ties and evolution of HMXBs in the SMC and in our Galaxy share 
a common nature and are governed by similar physical processes. 
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